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Abstract —A massive multiple-input multiple-output (MIMO) 
system, which utilizes a large number of base station (BS) 
antennas to serve a set of users, suffers from pilot contamination 
due to the inter-cell interference (ICI). In this letter, a graph 
coloring based pilot allocation (GC-PA) scheme is proposed 
to mitigate pilot contamination for multi-cell massive MIMO 
systems. Specifically, by exploiting the large-scale characteristics 
of fading channels, an interference graph is firstly constructed 
to describe the potential ICI relationship of all users. Then, with 
the limited pilot resource, the proposed GC-PA scheme aims to 
mitigate the potential ICI by efficiently allocating pilots among 
users in the interference graph. The performance gain of the 
proposed scheme is verified by simuiations. 

Index Terms —Massive muitiple-input multipie-output 

(MIMO), pilot contamination, graph coloring, pilot allocation. 

I. Introduction 

Massive multiple-input multiple-output (MIMO) has re¬ 
cently drawn extensive research interests due to its potential 
gains in spectrum efficiency and energy efficiency, whereby a 
base station (BS) equipped with a large number of antennas 
serves a set of users simultaneously 0, a- However, in 
practical multi-cell scenarios with unavoidable reuse of pilots 
in different cells, channel estimates obtained in a certain cell 
will be impaired by the reused pilots transmitted by users in 
other cells, which is referred as pilot contamination ITtt. It is 
widely recognized that pilot contamination is the performance 
bottleneck of massive MIMO systems ( 2 . 

How to mitigate pilot contamination has been extensively 
studied in the literature 0 - 0 . The greedy pilot allocation al¬ 
gorithm 0 can mitigate the pilot contamination by exploiting 
the statistical channel covariance information, but it suffers 
from high complexity. A blind method based on subspace 
partitioning 11 is able to reduce the inter-cell interference 
(ICI) when the channel vectors from different users are as¬ 
sumed to be orthogonal. An adaptive pilot clustering method 
0 based on coalitional game theory can mitigate the pilot 
contamination, whereby only a subset of the pilot resource is 
utilized in each cell, which leads to a slight loss in spectral 
efficiency. 

The graph coloring algorithms (GCAs) developed from the 
graph theory, which is a powerful mathematical tool, have been 
widely utilized to reduce ICI in mobile communications 0- 
0. For example, the available frequency bands are divided 
into several sub-bands and then allocated among femto-cells 
based on GCAs to strikes a balance between the ICI reduction 
and spatial frequency reuse 0, 0. A graph coloring based 
resource allocation scheme was proposed to reduce the ICI in 
the device-to-device (D 2 D) communication system by allocat¬ 
ing different resource to different D 2 D pairs with severe ICI 
0- These existing works try to mitigate the ICI by allocating 
different resource in the time or frequency domain to different 
femto-cells or D 2 D pairs based on GCAs, but all of them lead 
to a significant loss of spectral efficiency. 


In this letter, different from these existing works 0-0, 
which are essentially based on frequency division, we consider 
a multi-cell multi-user massive MIMO system where all cells 
work in the same frequency band to achieve higher spectral 
efficiency, and a graph coloring based pilot allocation (GC- 
PA) scheme is proposed to reduce the pilot contamination. 
Specifically, different from 0, 0 which construct the inter¬ 
ference graph among femto-cells, we construct the interference 
graph based the graphic structure of the ICI relationship among 
users in all cells, whereby all users in the same cell are 
connected to each other due to the intra-cell interference, and 
users in different cells having the potential ICI stronger than 
a certain threshold are also connected. After that, unlike 0, 
0 which simply allocate different resource to different femto- 
cells, the proposed GC-PA scheme greedily allocates different 
pilots to the connected users in the interference graph under 
the constraints of limited pilot resources and non pilot reuse 
within one cell. Simulation results verify the effectiveness and 
performance gain of the proposed GC-PA scheme in typical 
massive MIMO systems. 


II. System Model 

As shown in Fig. [j] (a), we consider a multi-cell massive 
MIMO system composed of L hexagonal cells, and each cell 
consists of a BS with M antennas and K (K <C M) single¬ 
antenna users 0 . The channel vector G C Mxl from 

the user (j, k) (k -th user in the j-th cell) to the )-th BS can 
be modeled as 


where denotes the large-scale fading coefficient, and 

g tj k),i ~ CM( 0 ,\m) denotes the small-scale fading vector. 

According to 0- 0, the pilots $ = [<p 1 <p 2 • • • 4>k\ T G 
C Kxt w j t j 1 length T (r > K) are orthogonal, i.e., = 

Ik, and the pilot resource is reused in adjacent cells. The 
classical pilot allocation scheme randomly allocates the pilot 
<fi to the user (j,k), and guarantees that every pilot will 
not be reused within one cell, i.e., p/j t k) 7^ P(j,k'),Vk 7^ k 
0. According to 0, the uplink signal-to-interference-plus- 
noise-ratio (SINR) of the user (j, k) can be calculated as 
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where aP k: , denotes the power of uncorrelated interference 
and noise which can be substantially reduced by increasing 
the number of BS antennas, V,., /??., ,,, . denotes 

pilot contamination caused by pilot reuse, and 1 ( :h k) denotes 
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Fig. 1. A typical example: (a) Classical random scheme: the available pilots will be randomly allocated to different users in each cell, which may cause severe 
ICI between users in adjacent cells, e.g., users in the shadow area; (b) Potential ICI: the ratio of interference channel strength and effective channel strength; 
(c) Construction of interference graph: an undirected graph can be constructed to describe the potential ICI relationship among all users; (d) Conventional 
GCAs: Totally C = 4 plots are required for LK = 9 users; (e) The proposed GC-PA scheme: only K = 3 pilots are utilized for LK = 9 users. 


the set of users with the same pilot as user (j, k), which is 
defined as = {(/,£/) : P(j>, k ') = P(j,k)} \ {(j,k)}- 

Thus, the corresponding average uplink achievable rate of the 
user ( j, k) can be calculated as 

C% k) = (1 - Mo)e{ log 2 (l + SINR^)}, (3) 

where po evaluates the loss of spectral efficiency caused by 
the uplink pilot transmission, which is actually the ratio of 
the pilot sequence’s length r and the channel coherence time 
l 0 , i.e., p 0 = 

It has been proved that the average uplink achievable rate 
is limited by pilot contamination (T). For example, as shown 
in Fig. [j] (a), the users in the shadow area suffer from severe 
pilot contamination due to the classical random pilot allocation 
scheme. It should be pointed out that above definitions and 
derivations can be directly applied in practical cellular network 
with irregular geometry (9). 

III. Graph Coloring based Pilot allocation 
A. Construction of Interference Graph 

From Eq. (0, it is clear that the ICI strength to a specific 
user (j, k) can be measured by E W,k')ei {j , k) Thus ’ 

we define the metric V{j,k),{j' ,fc') to measure the strength of 
potential ICI between user ( j, k) and ( j ', in different cells 
(j f j') when they are allocated with the same pilot, i.e., 

V{j,k),(j',k') = Plj-,k'),jlfi{j,k),j + $J,k),j'/P\j',k'),j'- W 

An intuitive illustration of V(;j.k),(j',k') * s shown in Fig. [T] (b), 
which indicates that rj is the ratio of the interference channel 
strength and the effective channel strength. Larger V(j,k),(j',k') 
indicates more severe ICI will be introduced between (j, k) 
and (j\ when the same pilot is allocated to them. 

In order to construct the interference graph to describe 
the potential ICI relationship among all users in the multi¬ 
cell massive MIMO system, a threshold 7* is introduced to 
determine whether two users in different cells can reuse the 
same pilot, i.e., V(j, k ),(j',k') > 7 th. where 7* will be discussed 
in detail in Section III-C. Based on 74 , we can generate a 
binary potential ICI matrix B = ,k')\LK xlk, where 

each entry EtUX/A') is defined as 

f 1 , j =/) k f k\ 

b U,k),U',k') = < 1. 3 > 7th, (5) 

[ 0 , otherwise. 

Intuitively, all users in the same cell are connected to each 
other, and two users in different cells are connected if their 
potential ICI strength is stronger than the threshold 7*. 


As illustrated in Fig. Q] (c), an interference graph can be 
constructed as an undirected graph G = (V, £) based on the 
potential ICI matrix B from Eq. (0. where the vertexes in 
set V denote all users, and the edges in set £ denote potential 
ICIs among users. 

B. Proposed GC-PA Scheme 

GCAs developed in the graph theory is able to color the 
vertexes of a graph with the minimum number of colors, C, 
such that no two connected vertexes will share the same color 
|| 6 l . Equivalently, for the pilot allocation problem in massive 
MIMO systems, GCAs can ensure that any connected users 
will be allocated with different pilots to mitigate the potential 
pilot contamination. Unfortunately, the minimum kinds of 
colors required by the conventional GCAs, i.e., C, is variable 
and unpredictable, since the interference graph consisted of 
moving users is dynamically changing due to the user mobility. 
For example, as shown in Fig. |T] (d), totally 4 kinds of colors 
are desirable to realize the graph coloring process, i.e., C = 4. 
However, in typical massive MIMO systems, only K = 3 
pilots are available to ensure orthogonal pilots used within one 
cell, and thus the conventional GCAs could not be applied. 


Algorithm 1 Proposed GC-PA Scheme 
Input: System parameters: K, L, and 7 th; 

Large-scale fading coefficients: 

Output: Pilot allocation: {pu.k) }■ 

1: iP(j,k)} = 0,0 = 0 , calculate 6 <j,k),</,&'>■ 

2 : while 3 P{j tk ) = 0 do 

3: (jo,ko) = argmax{ E b (j,k),(i',k')}- 

(j,k)£Q jy^j' 

4: A = {1, 2, • • • , I<} \ {p Uo , k) :l<k< I<}. 

5: A 0 = argmin { E b (joAo),(jA)}■ 

AeA PO,fc>= A 

6: P{ioM) = A 0 , O = O U (jo, k 0 ). 

7: end while 


In order to address this problem, the GC-PA scheme is 
proposed to mitigate the potential ICI with the constraint 
of predefined pilot resource for various user topology. The 
proposed GC-PA scheme greedily allocates the available pilot 
with smallest ICI to the most important user in a sequen¬ 
tial way, and the pseudocode is provided in Algorithm 1. 
Given the input information about cellular network parameters 
L. K, the threshold 7 ^, and the large-scale fading coefficients 
0(:j.k),ii this algorithm outputs the pilot allocation information 
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V(j.k) for a H users. The physical meaning of Algorithm 1 can 
be explained step by step as follows: 

1 ) Initialization (step 1 ): The pilot allocation is firstly 
initialized, i.e., {p( h k)} = 0 - Then, the set Ll containing 
the users which have been allocated with their corresponding 
pilots is also initialized as a null set, i.e., fi = 0. Given 
/, the binary potential ICI relationship t k),(j',k') can 
be calculated based on Eq. (0 

2 ) Loop Condition (step 2 ): The users will be allocated 
with their corresponding pilots in turns, and thus the loop 
condition is that there is at least one user who has not been 
allocated with any pilot, i.e., 3 P(j } k) = 0 . 

3 ) User Selection (step 3 ): Considering the fact that more 
connections in the interference graph lead to severer ICI, the 
user (jo , ko ), who has the maximum connections in adjacent 
cells and has not been allocated with any pilot, will be selected. 

4 ) Available Pilot Set (step 4 ): To avoid the intra-cell 
interference, no pilot will be reused within the same cell JT|, 
12 - The pilots, which have not been allocated to users in the 
jo-th cell, consists the available pilot set A for user (jo,ko). 

5 ) Pilot Selection (step 5 ): In order to allocate different 
pilots to connected users to mitigate ICI as much as possible, 
the pilot with index Ao, which has the minimal reused times 
among the connections of the selected user (jo,ko), will be 
selected from the available pilot set A for user (jo, kg). 

6 ) Pilot Allocation (step 6 ): In each loop, the selected pilot 
with index Ao will be allocated to the selected user (jo,ko), 
and then this user will be added into the set IT The procedure 
from step 2 ~ 6 will be carried out for KL times until all 
users are allocated with their corresponding pilots. 

In contrast to the conventional GCAs which requires an 
unpredictable number of pilots, the proposed GC-PA scheme 
is able to solve the pilot allocation problem with restricted 
K pilots. For instance as shown in Fig. [T] (e), total K = 3 
pilots are utilized to realize pilot allocation with significantly 
reduced ICI by the proposed GC-PA scheme. 

C. Analysis of Threshold y t h 

As the threshold 7* is used to be compared with the 
proposed heuristic metric it is a key parameter for 

the construction of the interference graph and the performance 
of the proposed GC-PA scheme. 

From Eq. ( 0 . we can find that any threshold chosen from 
( — 00,77min) leads to the same interference graph with the 
threshold 7* = r^min, where r/ m i n denotes the minimum value 
of 77 (j,k)Aj‘,k') based on Eq. 0 . Similarly, the upper bound 
of the threshold 7* can be set as rj max , so we have 

7th C [ 77 min j T^max] • (6) 

In one extreme case, when the threshold 7* is set to r) max , 
only users within the same cell will be connected to each other. 
In this case, both the GCAs and the proposed GC-PA scheme 
degrade to the classical random pilot allocation scheme, since 
no ICI information is exploited. By decreasing the threshold 
7th, more users in different cells will be connected to reflect the 
potential ICI relationship. In this case, as discussed above, the 
GCAs are able to eliminate the ICI among the connected users 


with an increasing number of pilots C, while the proposed GC- 
PA scheme can also mitigate most ICI with the restricted K 
pilots, where K < C. In the other extreme case, when the 
threshold 7* is set to ?7 m i n , all KL users will be connected to 
each other. It is evident that total C = KL pilots are required 
by the GCAs in this case, while the proposed GC-PA scheme 
degrades to the classical random pilot allocation scheme. 

To obtain a near-optimal threshold 7* in practical systems, 
we propose a simple yet effective method called iterative grid 
search (IGS) as 

7th = /igs ^ [ 77 min5 77 m ax] ■ N, t'J , ( 7 ) 

where N denotes the number of grids in each search step, and 
T denotes the number of iterations. Specifically, the interval 
[ 77 min, 77 max ] is uniformly sampled by N points in the first 
iteration, i.e., {7}^ : 1 < i < N}, 7^ = 77 min , 7^ = 
7 S -^ + A^ 1 ), and = r)m "~^ min . By denoting jml x as 

the one out of {7^} that can achieve the highest average 
uplink SINR, a sub interval [7max — A^/ 2 ,7max + A^ / 2 ] 
can be obtained after the first iteration, which will be further 
sampled in the next iteration to obtain a smaller range [7mlx — 
A^ 2 )/ 2 ,7mlx + A( 2 )/ 2 ]. This procedure is carried out in an 

iterative way for T times, and finally a near-optimal threshold 

(T) 

7th = 7max can be obtained. 

D. Discussion of Implementation 

In Section III-A, the interference graph is constructed ac¬ 
cording to the large-scale fading coefficients which 

actually denotes the channel strength between the users and 
the BSs. 

Compared with the existing works 0, 0 which assume 
the statistical channel covariance can be accurately acquired, 
the proposed GC-PA scheme only requires large-scale fading 
coefficients, which changes slowly and can be easily tracked 
with low complexity 0, 12, 0. For example, in practical 
mobile cellular networks of the long term evolution (FTE) sys¬ 
tems 0, an user firstly captures the cell-specific broadcasting 
signals to measure the channel conditions to available BSs, 
and then selects the BS with the best channel condition. After 
that, this user will continue to track the channel conditions 
to available BSs, and thus the handover among different cells 
can be easily realized. 

As addressed in FTE systems 0, the BSs in adjacent cells 
are connected to each other through X 2 interface to share and 
exchange information. Moreover, the coordinated BSs are all 
connected to a mobility management entity (MME) through 
S 1 interface, which has a powerful computing capability. Thus, 
the MME is able to collect the large-scale fading coefficients 
from coordinated BSs, and operate the proposed GC-PA 
scheme to provide a near-optimal pilot allocation pattern. 

According to Algorithm 1 , the complexity 0 (K 2 L 2 ) is 
required to operate the proposed GC-PA scheme based on 
a specific threshold 7*. According to the IGS procedure. 
Algorithm 1 will be carried out for NT times to achieve 
the near-optimal threshold 7*, and thus the total complexity 
is 0 (NTK 2 L 2 ), which is acceptable for the powerful MME 
0 . 
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Fig. 2. Simulation results: (a) The IGS procedure to obtain the near-optimal threshold 7 t h with L = 7, K = 8, M = 128, N = 20, and T = 2; (b) The 
CDF of the users’ uplink achievable rate with L = 4, K = 4, and M = 128; (c) The average uplink achievable rate per user against the number of BS 
antennas M with L = 7 and K = 8. 


TABLE I 

Simulation Parameters 


Number of cells L 

4, 7 

Number of BS antennas M 

10 < M < 10000 

Number of users in each cell K 

4, 8 

Cell radius R 

500 m 

Path loss exponent a 

3 

Log normal shadowing fading cr s h a dow 

8 dB 

Loss of spectral efficiency pq 

Mo = 0.2 

IGS parameters 

N = 20, T = 2 


IV. Numerical Results 

In this section, we consider a typical hexagonal cellular 
network with L cells, where each cell has K users with 
single-antenna and a BS with M antennas 12-El, and the 
default values of the various parameters are listed in Table 
U As addressed in JT], the large-scale fading coefficient 

P{j,k),i can be modeled as P(J,k),i = Z(j,k),i/(r(j,k),i/R) a , 

where z^ k),i denotes the shadow fading and possesses a log¬ 
normal distribution, i.e., 10 log 10 (^/ Jifc \ j) ~ CJV(0, er s hadow), 
and j is the distance from the user (j. k) to the i-th BS. 

Fig. [2] (a) shows the IGS procedure to obtain the near- 
optimal threshold 7 * when L = 7, K = 8, M = 128, and 
N = 20 are considered as an example, and we only select 
T as T = 2. The X-axis denotes the normalized threshold 
interval, and the Y-axis denotes the average uplink SINR of the 
proposed GC-PA scheme with the thresholds ( 7 ^ }, which are 
the uniformly sampled points in the i-th iteration. It is evident 
that the best sub-interval is selected after the first iteration of 
the IGS procedure, and a near-optimal threshold 7 * = 7 max 
can be obtained after the second iteration. 

Fig. U (b) plots the cumulative distribution function (CDF) 
curve of the users’ uplink achievable rate with L = 4, K = 4, 
and M = 128. The classical scheme randomly allocates the 
available pilots to all users JT), the conventional GCA is able to 
eliminate all potential ICI by consuming unrestricted number 
of pilots C 01, 0, while the optimal pilot allocation is 
obtained through exhaustive search among all possible cases 
as much as (Ff!) i_1 = (4!) 3 = 13824. The curve of the 
conventional GCA is denoted by a dotted line, since they 
could not be applied in practical systems with predefined and 
restricted pilot resource ID, H, ID. It should be pointed out 
that, pc = p 07 ^ instead of p 0 , is utilized to calculate the 
users’ uplink achievable rate of the conventional GCA. From 
Fig. |2](b), it is evident that the conventional GCA outperforms 
the classical scheme by about 0.6 bps/Hz, and another gain 
of about 0.4 bps/Hz can be achieved by the proposed GC-PA 


scheme. Moreover, the performance gap between the proposed 
GC-PA scheme and the optimal PA is only about 0.1 bps/Hz. 

Fig. [2] (c) presents the average uplink achievable rate per 
user against the number of BS antennas M with L = 7 
and K = 8. The curve of the optimal pilot allocation is not 
plotted here since its computational complexity is too high to 
be realized, i.e., (AT!) i_1 = (8!) 6 ss 4.3 x 10 27 . It is evident 
that the proposed GC-PA scheme outperforms the classical 
scheme by about 1.0 bps/Hz when the typical number of BS 
antennas M = 128 is considered. Moreover, it is noteworthy 
that the performance of the proposed GC-PA scheme can be 
continually improved by increasing the number of BS antennas 
M, while that of the classical scheme becomes saturated fast. 

V. Conclusions 

In this letter, we have proposed a pilot allocation scheme 
based on graph coloring to mitigate pilot contamination for 
multi-cell massive MIMO systems. Firstly, an interference 
graph is constructed to describe the potential ICI relationship 
among all users. After that, the GC-PA scheme is proposed to 
eliminate most potential ICI in the interference graph. 
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